Detection of Photometric Variations in the sdBV Star JL 166 
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ABSTRACT 

We report the discovery of oscillations in the hot subdwarf B star JL 166 from 
time-series photometry using the Goodman Spectrograph on the 4.1-m Southern 
Astrophysical Research Telescope. Previous spectroscopic and photometric ob- 
servations place the star near the hot end of the empirical sdB instability strip 
and imply the presence of a cool companion. Amplitude spectra of the stellar 
light curve reveal at least 10 independent pulsation modes with periods ranging 
from 97 to 178 s and amplitudes from 0.9 to 4 mma. We adopt atmospheric pa- 
rameters of Te// = 34350 K and log g = 5.75 from a model atmosphere analysis 
of our time-averaged, medium-resolution spectrum. 

Subject headings: stars: subdwarfs - stars: oscillations - stars: individual: JL 
166 



1. Introduction 

Hot subdwarf B (sdB) stars are a class of objects identified with models of extended 
horizontal branch stars and have temperatures ranging from 22000 to 40000 K and log g 
values from 5.0 to 6.0 dex. They dominate surveys of faint blue objects at high galactic 
latitudes and are often cited as the main source of the UV excess observed in giant elliptical 
galaxies. Although their origins and evolutionary tracks are still debated, the sdB stars are 
believed to be evolved, lower-mass (~O.5M0) stars with a He-burning core surrounded by a 
thin H layer. Models of sdB stars show they will evolve directly to the white dwarf cooling 



^Based on observations at the SOAR Telescope, a collaboration between CPNq-Brazil, NOAO, UNC, and 
MSU. 

^Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599-3255, USA; 
bbarlow@physics.unc.edu, bhdunlap@physics.unc.edu, clemens@physics.unc.edu 

■^Department of Physics, University of Oxford, Keble Road, Oxford 0X1 3RH, England; 
aelg@astro.ox.ac.uk 



- 2 - 



sequence after core He exhaustion. Their optical spectra are dominated by Balmer hues and 
sometimes also display He lines. 



Charpinet et al.l (Il996l . 119971 ) predicted the existence of a class of pulsating sdB stars 
based on the presence of p-modes in their model stars at the right temperatures and surface 
gravities. The oscillations they found were driven by a K-mechanism associate d with an opac 



i ty bu mp from the ionization of Fe. Contemporaneous with this prediction, Kilkenny et al. 



( I1997I ) reported the first detection of p-mode oscillations in the sdB star EC 14026-2647, 
opening up the possibility of using asteroseismological methods to probe the interiors of 
these stars. Since this discovery, more than 60 pulsating sdB stars (sdB Vs) have been ob- 
served with either p- or g-mode oscillations. The rapid p-mode pulsatorcl oscillate with 
periods between 80 and 600 s and amplitudes around 10 mma. The g-mode pulsator^ have 
periods ranging from 1 to 2 hrs and amplitudes comparable to the p-mode pulsators. In 
the log g-Teff plane, the p-mode oscillators cluster together in an instability strip with Tg// 
from 28000 to 35000 K and log g between 5.2 and 6.1 dex. The slow pulsators typically have 
cooler temperatures and lower surface gravities with Teff between 23000 and 30000 K and 
log g near 5.4 dex. In the region where the red edge of the p-mode instability strip overlaps 
the blue edge of the g-mode pulsator strip, three sdBVs have been discovered that show both 



p- and g-mode oscillations (ISchuh et al. 



20061 : iBaran et al.ll2006l : iLutz et all 120091 ). 



We report the discovery of p-mode oscilla t ions i n JL 166, a hot s ubdwarf B star that was 
first listed in the catalogue of lJaidee fc Lyngal (Il969l ). In their survey, IJaidee &: Lyngal (119691 ) 
used the Schmidt telescope of the Uppsala Southern Station at Mount Stromlo to search for 
faint violet stars in Southern galactic latitudes. 296 stars were catalogued in the survey; 
their colors were assigned "decidedly violet" or "possibly violet" designations, while their 
brightnesses were estimated as "bright," "intermediate," or "faint." JL 166 was classified 
as "decidedly violet" with "intermediate" brightness, and further observations with a 100 
cm reflector at Siding Spring Obse rvatory led to B and V measureme nts of 15.00 and 15.23 , 
respectively ( Jaidee fc Lynga|[l969 ). More recent 2MASS observations ( Skrutskie et al. I l2006h 
found J, H, and K magnitudes that imply the presence of a cool companion. The only 
published spectral data for JL 166 (IHebeii 119861 ) revealed strong Balmer lin es, He I lines , 
and the 4686 A He II line, and resulted in an sdOB spectral classification. iHeberl ( 19861 ) 
reported temperature and gravity values of 35500 K and 5.8 dex, respectively, which place 
JL 1 66 near the blue edge of the empirical sdB instability strip (see Fig. 1 of I Charpinet 
200lh . 



''referred to as V361 Hya or EC 14026 stars 

^also known as PG1716, V1093 Her, or Betsy stars 
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We first learned of JL 166 after searcliing tlie online Subdwarf Database fl0stensenll2OO6l ) 
for sdB stars in the instability strip. Noticing JL 166 fell in the p-mode instability strip, 
we observed the star to look for variations in the stellar luminosity. With the Goodman 
Spectrograph on the 4.1 m Southern Astrophysical Research (SOAR) telescope, we have 
detected at least 10 independent frequencies in the light curve with periods ranging from 97 
to 178 s and amplitudes from 0.9 to 4 mma. 



Heben (120091 ) reviews what is curre ntly kn o wn ab out sdB stars, a poorly understood late- 



stage of stellar evolution. In particular. iHeberl (120091 ) highlights the diversity and emphasizes 
(his §14) the need for more determinations of Teff and log g. Further studies, such as the 
one on JL 166 presented here, would in due course permit a better, informed study of sdB 
star evolution and their contribution to the ultraviolet upturn seen in the spectra of giant 
elliptical galaxies. 



2. Photometry 

We obtained time-series photometry for JL 166 with the Goodman Spectrograph on the 
4.1-m SOAR telescope. In imaging mode, the camera-collimator combination re- images the 
SOAR telescope focal plane with a focal reduction of three times, resulting in a plate scale 
of 0.15 arcsec pixel"^. The camera houses a 4k x 4k Fairchild back-illuminated CCD with 
electronics and dewar provided by Spectral Instruments, Inc. Using optics of fused silica, 
Na Cl, and CaF2, the ent ire system is optimized for high throughput from 320 to 850 nm. 



See IClemens et al.l (120041 ) for further details on the spectrograph. 



We observed JL 166 on two engineering nights in 2008 September and October, obtaining 
uninterrupted runs of over 2 hours on each night. Each run was obtained through a red- 
blocking S8612 filter, which has a bandpass of 300 to 700 nm. To minimize the processing 
time between exposures yet keep the readnoise low, we read out a small subsection of the 
CCD (approximately 500 x 500 pixels binned 2x2) and used a readout speed of 100 kHz. 
Each exposure had an integration time of 10 s, resulting in a cycle time of ~13 s. Table 1 
summarizes the details of our photometric observations. 

In order to obtain light curves from the raw images, we extracted our photometry using 
the external IRAF package CCD_ HSP developed by Anton io Kanaan, which employs the 



aperture photometry preferred by lO'Donoghue et al.l (120001 ). Aperture widths were chosen 



to maximize the signal-to-noise ratio in the light curves and were approximately 1.7 times 
the seeing width. To subtract the sky from each stellar aperture, we used sky annuli with 
widths of 1.5 arcsec that started approximately 4.5 arcsec from the centers of the stars. 
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We divided our light curves by those of a constant comparison star to remove small-scale 
variations in the sky transparency. The comparison star used is located approximately 02' 
39'.'4 southwest of JL 166 and was the only one present in our frames with an adequate signal- 
to-noise ratio. Atmospheric extinction effects were corrected by fitting and normalizing the 
curves with parabolas. We note that this normalization may remove real variations in the 
stellar brightness on the order of our run length. Due to superb telescope guiding, we did 
not flat-field or bias-subtract any of the frames. The reduced light curves for JL 166 are 
presented in Figure 1. 

We analyzed our reduced light curves by combining Fouri e r ana lysis and least-squares 
fits in a standard manner using Period04 (see iLenz fc Bregerl (120051 ) for program details). 
As our two observing runs were separated by a month, the light curves were analyzed on a 
night-by-night basis. The amplitude spectra are displayed in Figure 2. We used a prewhiten- 
ing technique to perform our temporal analysis, fitting the largest signal in the amplitude 
spectrum, subtracting the fit from the data, and re-calculating the Fourier transform of the 
residual light curve. We repeated this process until the ca ndidate peak s had confidence levels 
below 99%, as given by the statistical test proposed by iKoeru (1l990l ). After each iteration 
we applied a non-linear, least-squares fitting routine to simultaneously fit the periods, am- 
plitudes, and phases of the isolated frequencies. Figure 3 shows the amplitude spectra at 
various stages of this iterative process. The extracted frequencies are printed in Table 2 and 
are each labeled fj, where j is ordered in terms of decreasing average amplitude. The errors 
shown are derived from the least-squares fits ; the actual errors may be larger by a factor of 



three (see [Montgomery &: O'Donoghud Il999l ) 



The multiperiodic nature of JL 166 is apparent from the amplitude spectrum of the 
complete light curve. We report 10 independent oscillations with periods ranging from 97 s 
to 178 s and amplitudes from 0.9 to 4 mma; the least significant of these has a false alarm 
probability near 10~^. Some of the reported frequencies were only detected on one of the 
two nights. Amplitude variability is dete cted in many sdB stars observed over an extended 
period; see for example iReed et al.l (120071 ). Whether the non-detection of some of the modes 
on one of the two nights is due to the amplitudes decreasing below detection limits is unclear 
since the amplitudes of the majority of the modes varied between the observation sets. Most 
of the power in the spectrum is concentrated between 6800 and 7600 //Hz. The highest- 
frequency variation (13461 ^Hz) is a combination frequency of the dominant peak, /i, and 
another lower- amplitude frequency, fg. Interestingly, neither fg nor the cross-frequency is 
present in the October data. 



We note that additional oscillation modes may be present in the amplitude spectrum. 
The mean noise level in the 4 to 14 mHz range is greater than that on either side of this 
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band. This effective plateau in the amphtude spectrum is visible both before and after 
prewhitening and suggests the presence of additional frequencies that cannot be identified 
in our data. 



3. Spectroscopy 



To improve upon the spectroscopy of lHeberl (jl986l ). we obtained four medium-resolution 
spectra of JL 166 with a combined integration time of 1800 s using the Goodman Spectro- 
graph on 2008 September 19. We used a 1.03 arcsecond slit and the 600 mm~^ VPH grating 
(0.65 A pixel"^ dispersion) to cover a spectral range from 3500 to 6300 A with a resolution 
of 4.5 A. Our spectral images were binned 2x2, resulting in 3.4 pixels per resolution element 
and a signal-to-noise ratio of approximately 70, measured at the continuum near 4200 A. 
Table 3 summarizes the details of our spectroscopic observations. 

We reduced our spectra with IRAF using standard bias-subtraction and wavelength- 
calibration procedures. A standard flat-field correction could not be applied due to a lack of 
proper fiat frames. We corrected this potential issue by applying a flux calibration using a 
spectrophotometric standard star (EG 21) whose spectrum was aligned on the same pixels as 
that of the JL 166 spectra. The reduced spectral frames were averaged to produce the single, 
time-averaged spectrum for JL 166 shown in Figure 4. As expected, our mean spectrum is 
dominated by Balmer li nes, He I line s, and the 4686 A He II line, confirming the sdOB 
classification assigned by iHeberl (119861 ). We deriv e a red shift of 79 ± 27 km s~^ from these 
lines, also in agreement with the results of iHeben (119861 ). 



Although typical indications of a low- mass main sequence companion (Ca I triplet, Mg 
I triplet, G band) are absent from the spectrum, Na D absorption lines (5890 A and 5896 
A) are present and confirm the existence of a cool companion (see inset of Figure 4). An 
over-subtraction of the Na D sky emission features is unlikely as the removal of the 01 sky 
emission line at 5577 A left no residuals. Moreover, the Na features in our spectra have 
redshifts comparable to those computed from the H and He features, and, consequently, we 
believe they are real. Contamination from interstellar Na absorption cannot be completely 
ruled out but is unlikely due to the high galactic latitude (-69.8°) of JL 166. As no other 
spectral features typical of a cool companion were observed, we canno t classify the companion 
using our data. The optical-IR color-color plots of I St ark et al.l (12004J ) suggest the companion 
is an M-type star contributing less than 5% of the total flux (in the V-band) from the system. 

To determine the atmospheric pa rameters of JL 16 6, we employed the model grid of 
synthetic optical spectra described in iHan et al.l (120071 ). The grid is defined in terms of 
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21 values for the effective temperature from 20000 K to 40000 K (in steps of 1000 K), 11 
values for the surface gravity from log g of 5.0 to 7.0 dex (in steps of 0.2 dex), and 4 values 
of the helium-to-hydrogen ratio from log N{B.e)/N{}l) of -3 to dex (in steps of 1 dex). 
Further resolution in the models was achieved through a standard grid interpolation routine 



(jPress et al.l 120071 ). We degraded each model spectrum by convolving it with a Gaussian 
whose FWHM mimicked that of a single resolution element in our spectrum. After re- 
binning the convolved models (while conserving flux) to match the binning of our data, we 
fit the available absorption lines in our mean spectrum using a least-squares analysis. 

No parameter set consistent with both the H Balmer lines and He ionization equilibrium 
could be obtained. This so-called "Balmer line problem" is not uncommon amongst sdB stars 
(iHeber et al.lll999l ). The simultaneous fitting of H Balmer lines results in Teff = 33300 K and 



log g = 5.75. Since the Hel/Hell ratio is highly sensitive to temperature and the He lines are 
relatively unaffected by gravity, we can also determine Teff by simultaneously fitting the He 
lines from both ionizations states at the gravity fixed by the Balmer line fits. This ionization 
equilbrium indicates Teff = 35400 K and log A/^(He)/iV(H) = -0.8. As a compromise, we 
adopt the average parameters deduced from the He ionization equilibrium and the H Balmer 
lines: Te// = 34350 ± 1000 K, log g = 5.75 ± 0.23, and log N{Re)/N(}l) = -0.8 ± 0.3. 
Figure 5 shows the best fits to the H and He profiles. 

To illustrate the discrepancies in the model fits. Figure 6 presents the gravity and 
temperature values resulting from the individual fits to the H Balmer lines, the simultaneous 
fit to the H Balmer lines, and the He ionization equilbrium requirement. The H Balmer line 
fits span a temperature range from 32300 to 33200 K and a log g range from 5.4 to 6.3 dex. 
A trend of increasing log g with increasing excitation level is apparent in the model fits. 
The temperature fits, on the other hand, display no obvious trend. The simultaneous fit to 
the H/? to Hr] lines shows that the temperature fit is most regulated by the lower-excitation 
lines. The location of the gravity fit is near the average value but is slightly partial to the 
higher-excitation lines. As previously noted, the He ionization equilibrium condition requires 
a temperature around 2000 K hot ter than that resulting from the Balmer line fits. This result 



is consistent with the findings of iHeber et al.l (119991 ) and has yet to be clearly resolved 



4. Discussion 

We have discovered a new member of the sdBV pulsator class, JL 166, from high signal- 
to-noise photometry. Our spectroscopic observations confirm a sdOB classification and find 
the spectrum best fits a hot subdwarf model with Tg// = 34350 ± 1000 K, log g = 5.75 ± 0.23 
dex, and log N{}ie)/N{}i) = -0.8 ± 0.3. These parameters are in agreement with those cited 
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by iHeberl (119861 ) . The infrared excess in the J, H, and K bands and our detection of weak 
Na D hnes suggest the presence of a cool companion, probably an M-type star. Although 
the contamination from the companion appears to be low (Mg I, Ca I, and the G-band are 
not visible), we note that our spectroscopic fits are limited by imcomplete knowledge of 
contamination by the companion. 

The analyses of our photometric observations show that JL 166 is a multimode pulsator, 
with at least 10 independent frequencies detected between and 5601 /iHz (178 s) and 10344 
/iHz (97 s) with amplitudes from 0.9 to 4 mma; the most dominant signal is at 134 s. These 
oscillations are consistent w ith p-mode pulsat ions in hot subdwarf B stars and appear to 
follow the log ^-P relation of iKoen et al.l (119991 ). Moreover, the frequencies, amplitudes, and 
number of the observed modes seem comparable to other sdBVs with similar atmospheric 
parameters. PG 1047+003, most notably, ha s Teff and log g valu es comparable to JL 166 
and pulsates with periods from 104 to 175 s (IKilkenny et al.ll2002l ). 



The amphtudes of the detected frequencies differ between observation sets. These varia- 
tions are either intrinsic to the star or are due to beating effects from unresolved frequencies. 
The residuals in the amplitude spectra after prewhitening suggest the presence of additional, 
unresolved frequencies. Intrinsic amplitude variations cannot be ruled out either, though, 
as they are observed in other sdBVs. Some of the oscillation modes in PG 1047+003, 
for example, undergo co nsiderable amplitude variations over a timescale of a few weeks 
(jO'Donoghue et al.lll998l ). More intensive photometric campaigns are needed to increase the 
resolution in the amplitude spectrum and distinguish between the two possibilities. 

Cross-correlation of the four spectra obtained showed the radial velocity of the JL 166 
sdB component to have changed by 0.0 ± 5 km s~^ in 0.02 days. No spectroscopic evi- 
dence which would justify interpretting the /i frequency change, between 2008 September 
17 and October 14, as a Doppler shift is therefore presently available. Furthermore, the 
corresponding change in /2 could not be explained by a Doppler shift which might explain 
the /i variation. 
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Table 1: Time-series Photometry Log 



UT Date 


Start Time 


T 

-■- exp 


T 

-■- cycle 


Length 


Filter 


(2008) 


(UTC) 


(«) 




(hr) 




Sep 17 


00:57:56 


10 


13 


2.21 


S8612 


Oct 14 


00:00:00 


10 


13 


2.22 


S8612 
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Table 2: Periodicites detected in the light curve of JL 166. 



ID 


UT Date 


Period 


Frequency 


Amplitude 




(2008) 


(s) 


(/iHz) 


(mma) 


/I 
/I 


Sep 


i / 


io4.zi lb U.Uo 


/^4ol.U ± 6.1 


4.Uo ± U.i / 




Uct 


1 A 


133.41 ± 0.09 


7495.7 ± 4.9 


3.33 ± 0.22 


/2 


Sep 


1 ^ 
i / 


146.40 ± 0.34 


6831 ± 16 


c\ r\r\ 1 r\ 

2.90 ± 1.0 




UCt 


1 A 

14 


144.08 ± 0.31 


6941 ± 15 


1.58 ± 0.22 


/3 


bep 


1 T 


131.36 ± 0.16 


7612.8 ± 9.5 


1.35 ± 0.17 




Oct 


14 


130.66 ± 0.13 


7653.5 ± 7.6 


2.12 ± 0.22 


/4 


Sep 


17 


140.82 ± 0.14 


^ -x r\'\ 'X 1 ^7 f~\ 

7101.1 ± 7.0 


1.85 ± 0.18 




Oct 


14 


141.07 ± 0.29 


7089 ± 14 


1.49 ± 0.23 


/5 


Sep 


17 


147.84 ± 0.52 


6764 ± 24 


1.89 ± 1.0 




Oct 


14 


147.30 ± 0.33 


6789 ± 15 


1.42 ± 0.22 


/6 


Sep 


17 










Oct 


14 


96.68 ± 0.10 


10344 ± 10 


1.48 ± 0.21 


/7 


Sep 


17 


101.42 ± 0.12 


9860 ± 12 


0.99 ± 0.17 




Oct 


14 


101.38 ± 0.11 


9864 ± 10 


1.46 ± 0.21 


/8 


Sep 


17 










Oct 


14 


178.53 ± 0.40 


5601 ± 13 


1.20 ± 0.21 


/9 


Sep 


17 


167.35 ± 0.24 


5975.5 ± 8.5 


1.06 ± 0.17 




Oct 


14 








/lO 


Sep 


17 


109.52 ± 0.16 


9131 ± 13 


0.91 ± 0.17 




Oct 


14 








/I +/9 


Sep 


17 


74.48 ± 0.05 


13426.6 ± 9.0 


0.86 ± 0.17 




Oct 


14 
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Table 3: Spectroscopy Log 



Date 


T 


Grating 


Slit 


Resolution 


(HJD) 


(s) 


(1/mm) 


(arcscc) 


(A) 


2454728.75552 


300 


600 


1.03 


4.5 


2454728.76029 


500 


600 


1.03 


4.5 


2454728.76708 


500 


600 


1.03 


4.5 


2454728.77381 


500 


600 


1.03 


4.5 



500 



1500 



2000 



0.015 


-0.015 
0.015 


-0.015 
0.015 



0.015 


-0.015 









0.015 













-0.015 






' ' ' ' 1 ' ' ' ' 1 ' ' ' 1 


1 1 ■ 1 ' : 


Amplitude 


0.015 


i i i i 1 1 i i 


i 1 i i i i 1 1 i i 1 : 




. * **^ 




-0.015 






1 1 1 1 1 1 i 1 1 1 1 1 


I'll 

• •% .V '*. 


Fractional 


0.015 


-0.015 


1 1 1 1 1 ' 1 1 
,•* - ' •• . 


1 1 i 1 1 1 1 1 1 1 1 : 

■•• - ■' ? ■ .. " 


1 1 1 i 1 ' ' i i 1 i 1 1 1 


I'll' 




0.015 



1 1 i i 1 1 1 i 


1 1 1 1. 1 i 1 1 1 1 1 : 
■ -'' • •• •' ■. ■. •. 


. . . . 1 .... 1 ... . 


. . . . 




-0.015 


1 . . . 


1 .... 1 ... . 



6500 7000 7500 

Time from Run Start (s) 



6000 



6500 7000 7500 

Time from Run Start (s) 



8000 



Fig. 1. — Light curves for JL 166 from 2008 September 17 (left panel) and October 14 
(right panel). Small-scale variations in the sky transparency were removed by dividing the 
light curves by those of constant comparison stars. The light curves were then normalized 
with parabolas in order to remove atmospheric extinction effects. The data are presented 
unsmoothed. 
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Fig. 2. — Fourier amplitude spectra for the complete light curves from 2008 September 17 
(top) and October 14 (bottom) in the 0-35 mHz bandpass. 
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Fig. 3. — Amplitude spectra from 2008 September 17 (left panel) and October 14 (right 
panel). Spectra are shown before prewhitening (top), after subtracting the 4 largest peaks 
(middle) , and after subtracting the 9 largest peaks (bottom) . The subtracted peaks in each 
spectrum are marked in the preceeding amplitude spectra. Only signals with confidence 
levels above 99% (approximately indicated by the dashed lines) are claimed. 
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Fig. 4. — Time-averaged, flux-calibrated spectrogram of JL 166. The spectrum is an average 
of four 60-s exposures and was flux-calibrated using the spcctrophotometric standard star 
EG 21. The inset shows a magniflcation of the reddest portion of the spectrum in which the 
Na I (D) hues at 5890 A and 5896 A are visible next to a He I line as a single, unresolved 
feature. 
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Fig. 5. — Model fits (dotted lines) to the H Balmer and He line profiles in the time-averaged 
spectrum (solid hnes) of JL 166. 
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Fig. 6. — Positions of the best individual fits to the H Balmer hnes (denoted by their 
hne designations), the best simultaneous fit to the H Balmer lines (open circle), and the 
He ionization equilibrium solution (filled square) in the log g-Teff plane. Our compromise 
solution (cross) is also shown and has Tg// = 34350 K and log g = 5.75. 
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